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' We compute the black hole horizon entanglement entropy Se for a 

massless scalar field, first with a hard cutoff and then with high fre- 
quency dispersion, both imposed in a frame that falls freely across the 
horizon. Using WKB methods, we find that Se is finite for a hard cut- 
off or super-luminal dispersion, because the mode oscillations do not 
diverge at the horizon and the contribution of high transverse momenta 
' is cut off by the angular momentum barrier. For sub-luminal dispersion 

the entropy depends on the behavior at arbitrarily high transverse mo- 
menta. In all cases it scales with the horizon area. For the hard cutoff 
. it is linear in the cutoff, rather than quadratic. This discrepancy from 

the familiar result arises from the difference between the free-fall frame 
and the static frame in which a cutoff is usually imposed. In the super- 
luminal case the entropy scales with a fractional power of the cutoff 
that depends on the index of the dispersion relation. Implications for 
the possible relation between regularized entanglement entropy and the 
Bekenstein-Hawking entropy are discussed. An appendix provides an 
explicit derivation of the entangled, thermal nature of the near- horizon 
free fall vacuum for a dispersive scalar field in four dimensions. 
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1 Introduction 



All indications suggest that the Bekenstein-Hawking entropy of a black hole 



is true thermodynamic entropy [TJ [2J [3]. (Here A is the area of the horizon 
and L p the Planck length.) At the microscopic level it must therefore corre- 
spond to the logarithm of the number of thermodynamically relevant states. 
One can write down a formal path integral expression for the partition func- 
tion and entropy of the quantum gravitational statistical ensemble[4j, which 
in the classical approximation yields precisely this entropy, however the mi- 
croscopic nature of the states thereby counted remains obscure. 

It would appear that an unavoidable contribution to the entropy budget 
is the entanglement entropy Se across the horizon of quantum fields propa- 
gating on the black hole background [5J [6]. The entanglement entropy arises 
from the spacelike correlations in the freely falling vacuum (giving rise to 
Hawking radiation), on account of which the restriction to the region out- 
side the horizon is a mixed state. This can also be thought of as entropy of 
the thermal atmosphere of the black hole [ZJ |8], particularly referring to the 
near-horizon region where most of the entropy is localized. It is difficult to 
see any rationale for omitting the entanglement/atmosphere entropy. How- 
ever, in the standard semiclassical evaluation it diverges due to the infinite 
density of states at the horizon. 

One can only conclude that the divergence in Se must be somehow cut 
off by quantum gravity effects. A simple dimensional analysis of the geom- 
etry suggests that when cut off at Planck scale, the entanglement entropy 
across the horizon will scale as A/L 2 , just like the Bekenstein-Hawking en- 
tropy. If so, then Se would make an important contribution to Sbh- This 
presents a "species problem," since Se depends on the number and nature 
of quantum fields, whereas the Bekenstein-Hawking entropy is universally 
given by A/4L 2 ,. 

A natural solution to this problem may exist however p2 [10]: the value 
of the Planck length L p = (/iG/c 3 ) 1 / 2 appearing in Sbh involves the low 
energy gravitational constant G, whose value should depend, though renor- 
malization effects, on the quantum matter fields. The shift of 1/G produced 
by the matter fields could induce a shift of Sbh that corresponds precisely 
to Se- The formal arguments of Refs. [10] support this possibility, as do 
explicit regularized computations for free field theory such as in [11]. How- 
ever, mismatches seem to arise for contributions of vector fields [12] and 





2 



in odd dimensional spacetimes [13]. On the other hand, very general argu- 
ments |144 [To] indicate that if one focuses not on entanglement entropy but 
only on the contribution of matter to the entropy derived from the thermal 
partition function for gravity plus matter (which may be the same thing in 
some cases), then the relation between black hole entropy and renormaliza- 
tion of G is universal. A purely thermodynamic argument [16] suggests the 
even stronger statement that G is wholly determined by the entropy of local 
Rindler horizons. 

Previous studies of black hole entanglement entropy with a physically 
motivated regulator (as opposed to a technical regulator such as Pauli- 
Villars) have adopted a cutoff scheme where a highest momentum [6] or 
shortest distance [7] was imposed in the static frame. In both cases it was 
found that the entropy scales as A/e 2 where e is the cutoff length. However, 
the equivalence principle suggests that if there is a cutoff on the freely falling 
vacuum imposed by quantum gravity it should be physically determined in 
a freely falling frame. 

The aim of this paper is to investigate the nature of the regulated entan- 
glement entropy when a freely falling momentum cutoff A is imposed. We 
shall consider both a hard momentum cutoff, and a cutoff imposed by the 
introduction of dispersion at high wavevectors, and we carry out the analysis 
using WKB methods. We find that unlike what is found in the static frame, 
the entropy scales linearly with the cutoff as An v A for a hard cutoff, and 
with a similar but more complicated scaling law in the case of dispersion. 
In this expression k v is the surface gravity measured in the freely falling 
frame in which A has been defined. When that frame is at rest at infinity, 
k v reduces to the standard surface gravity. 

A surprising conclusion may be drawn from our results. If the freely 
falling cutoff A is identified with the Planck mass, then the regulated en- 
tanglement entropy is negligibly small compared to the Bekenstein-Hawking 
entropy for macroscopic black holes. This would provide a different resolu- 
tion to the species problem mentioned above (unless the number of species 
were of order M p /k v ). On the other hand, as argued in |17[ 118] . it is also 
possible that the cutoff arises from quantum gravitational effects in the black 
hole background, and that its value equals the Planck mass when measured 
in the static framell Such a cutoff would appear as A ~ M 2 /k v in the free- 
fall frame, and the regulated entropy would be of order A/L 2 after all. We 

3 If quantum horizon fluctuations fuzz out the horizon over the much longer scale 
(RsLp) 1 / 3 , as suggested in [191 1201 I21j . then the question arises whether the entropy 
is cut off at this scale or instead if the cutoff "floats" with the fluctuations. 
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shall return to a discussion of this issue at the end of the paper. 

Whatever the regularization scheme one chooses, the divergence in the 
entanglement entropy can be traced to two distinct sources. First, the den- 
sity of modes at fixed frequency and angular momentum diverges due to 
an infinite number of oscillations as the horizon is approached. Second, 
there is the sum over all angular momenta. In our approach, we find that 
in the case of both the hard cutoff and dispersion the WKB mode density 
is regulated. In effect, the dispersion introduces a maximum radial wave 
number contributing to the thermal atmosphere of the black hole. We fur- 
ther find that the mode density vanishes at a finite angular momentum in 
the case of a hard cutoff or super-luminal dispersion, while it decreases to 
zero asymptotically at large angular momenta in the sub-luminal case. The 
angular momentum sum is therefore finite in the former cases, while in the 
latter case it depends on the dispersion at asymptotically large momenta 
and cannot be determined within the WKB approximation because of the 
asymptotically vanishing mode density. 

The effect of sub-luminal dispersion on entanglement entropy has been 
previously investigated in Ref . [22] . There the dispersion was introduced in 
the static (accelerated) frame of the black hole, rather than in the freely 
falling frame. In that case, there is a piling up of transplanckian modes at 
some radius greater than the horizon radius, so the mode density at fixed 
angular momentum remains infinite and therefore the entropy diverges. This 
is not the result reported in [22J, where it was instead concluded that the 
entropy is finite. We believe this conclusion was based on an incorrect 
evaluation of the entropy. 

In following sections of this paper we first review the nature of the diver- 
gent entanglement entropy in the usual relativistic case. Next we introduce 
a hard momentum cut-off in the freely falling frame and analyze how it reg- 
ulates the divergences. Finally we introduce dispersion at high momenta 
also in the freely falling frame, and analyze its impact on the divergences 
of the entanglement entropy. A long appendix provides details on the field 
modes, the entangled, thermal nature of the free fall vacuum, and the WKB 
evaluation of the mode density governing the entanglement entropy. 

2 Thermal atmosphere of a black hole 

We consider nonrotating, static black hole spacetimes and adopt Painleve- 
Gullstrand (PG) coordinates since they are adapted to the free-fall frame 
in which the momentum cutoff or dispersion will be introduced. A free-fall 
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frame is defined by a family of radial infalling geodesies all with the same 
velocity at infinity. The corresponding PG line element in four dimensions 
takes the form 

ds 2 = dr 2 - (dp - v(p)dT) 2 - r(p) 2 (d6 2 + sin 2 Odcf 2 ), (2) 

in units with the speed of light c = 1 (see e.g. Ref. [23]). The spacetime has 
r-translation and rotational symmetry. The coordinate r differs from the 
Schwarzschild-like time by a constant scaling and an additive function of p. 
The function v(p) is the velocity dp/dr of the radial geodesies, along which 
r measures the proper time. These geodesies are orthogonal to the surfaces 
of constant r, on which p measures radial proper distance. We consider the 
case v(p) < 0, so that these geodesies are infalling. A surface p = pn on 
which v 2 (ph) = 1 is lightlike and is a future event horizon, i.e. no timelike 
curve can exit the interior region p < pjj- 
Expanding around the horizon one has 

v(p) = -1 + n v (p - p H ) H . (3) 

As we shall see, the non-linear terms are not relevant to the entropy. The 
constant k v depends on the free-fall frame at the horizonQ It is equal to the 
surface gravity defined with respect to the Killing vector <9 T , which coincides 
with the standard surface gravity k of the horizon in the case that v(p) — > 
as p — > oo. For spacetimes in which v(p) approaches a nonzero constant Voo, 
d T has norm (1 — v 2 ^) 1 ^ 2 at spatial infinity, so 

^ = K (1-0 1/2 - (4) 

For the example of a Schwarzschild black hole of mass M, the frame that 
falls from rest at infinity has v(p) = — (2GM/p) 1 / 2 . In this case the surface 
gravity is k = and the area radius r(p) is equal to p, the radial 

proper distance measured at constant rH 

When a black hole forms from collapse of matter, the quantum state of 
any relativistic field tends towards the Unruh vacuum [2H [25] . This sta- 
tionary state coincides with the vacuum as defined by freely falling (FF) 

4 A coordinate independent local expression for k v is (l/2)(s ■ V£ 2 )/(s ■ £) 2 , where s is 
the outward radial unit spacelike vector orthogonal to the free-fall worldline at the horizon 
and £ is the Killing vector with arbitrary normalization (which cancels out in this ratio). 

5 The vanishing of r"(p) is equivalent to the vanishing of the (ingoing or outgoing) 
radial null- null component of the Ricci tensor. The condition r'(p) = 1 holds only if v(p) 
vanishes at infinity. 
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observers for modes with FF frequencies much larger than the surface grav- 
ity k. In this state fluctuations on either side of the horizon are corre- 
lated in such a way that when tracing over the inside degrees of freedom, 
the reduced density matrix for the outgoing modes has the thermal form 
p cxt = Z^ 1 exp(—H/Tjj), where H is the Hamiltonian generating the Killing 
time translations (normalized at infinity) and Tjj = k/2it is the Hawking 
temperature [261 EH EE]- (Here and below we adopt units with h = c = 1.) 
Some of the outgoing thermal quanta reach spatial infinity and constitute 
Hawking radiation. Most of them however fall back across the horizon into 
the black hole, forming a thermal "atmosphere" [29J. The entropy of this 
atmosphere S a tmo consists of the part of the entanglement (von Neumann) 
entropy — Ti(p ext In p e xt) contributed by the near-horizon modes, that is, not 
including the entropy of the few low angular momentum Hawking quanta 
that escape to infinity. Satmo is in fact dominated by initially outgoing 
modes with large angular momenta which are turned back by the angular 
momentum barrier very close to the horizon. 

Given the symmetries and the thermal character of the reduced density 
matrix p cxt , the atmosphere entropy can be expressed as a mode sum in the 
form 

/oo 
dw^2(2t + l)n(u,£)stheTm(u/T H ), (5) 

1=0 

where u> is the Killing frequency, t is the angular momentum, Stherm(kV^if) 
is the thermal entropy per modejl and n(u>,£) is the density of outgoing 
modes entangled to inside configurations. 

The entropy diverges for two independent reasons. First, for each fixed 
uj and £, the mode density has a UV divergence logarithmic in the radial 
momentum cut-off. This arises from the infinite number of mode oscillations 
caused by the infinite blue shift at the horizon. Secondly, the sum over 
angular momenta is unbounded, producing a UV divergence quadratic in 
the angular momentum. Notice that for a 1+1 dimensional black hole only 
the first, logarithmic divergence is present. 

6 Strictly speaking there is a weak ^-dependence of the mode entropy due to the trans- 
mitted part of the infalling modes. These start from infinity in the ground state and 
partially penetrate the effective potential barrier, whereas the reflected parts of the near- 
horizon outgoing modes are thermally populated. Since the transmitted part vanishes for 
I S> uj/k, this I dependence can be safely neglected when computing the entropy of the 
atmosphere. 
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3 Freely falling cut-off 



We are interested in regulating the atmosphere entropy by modifying the 
field propagation at high momenta as defined in a freely falling frame. This is 
in contrast to 't Hooft's "brick wall" approach |7J to regulating the entropy, 
where the cut-off is defined on static time slices. We first consider the simple 
case of a hard momentum cut-off. 

Since the momentum cut-off is stationary and spherically symmetric, 
the atmosphere entropy still has the form ([5]). The entropy per mode is 
unchanged, but the mode density n(u, £) is cut-off. To compute the atmo- 
sphere entropy we therefore only need to evaluate the cut-off mode density 
n cut( w > The calculation is identical for any free-fall frame in which the 
cutoff is imposed, so for notational simplicity in the following we adopt the 
one for which Voo = 0. 

We consider a massless scalar field in the spacetime with line element 
([2]). The field theoretic justification for our evaluation of the entropy, while 
essential for a full understanding, is sequestered in the Appendix so as not 
to obscure the conceptual structure and simplicity of the entropy evaluation 
via WKB methods. The reader may wish to read the Appendix at this point 
before coming back to the main flow of the entropy computations. 

The modes of fixed frequency and angular momentum take the form 

cputm = exp(-iu^)r -1 u w< (p)li m (0, <f>). (6) 

Since the relevant modes have very short wavelength, the radial mode func- 
tions u U} £(p) are accurately described by the WKB form u ~ exp(i / kdp), 
where the radial momentum k = k 0J i(p) obeys the relativistic dispersion 
relation g^k^kp = 0, which here takes the form 

(uj-v(p)k) 2 = k 2 +£(£ + l)/r 2 . (7) 

The quantity f2 = lo— v k is the free-fall frequency, whereas u is the conserved 
Killing frequency. We introduce the UV cutoff by counting only modes 
with k < A for some fixed A 3> k. The cutoff is imposed on the radial 
component of momentum only, while the transverse (angular) momentum is 
unrestricted. It turns out however that the result would differ little if the 
cutoff were imposed on the total momentum, since the angular momentum 
barrier imposes an effective transverse momentum cut-off of order (wA) 1 / 2 <C 
A. 

The atmosphere entropy is defined as the entropy of the density matrix 
obtained after tracing over the field degrees of freedom localized inside the 
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horizon. We shall evaluate the entropy here with the help of the WKB 
approximation, in which the modes are characterized by a well-defined mo- 
mentum k^e^p), the solution of (|7|), at any radial position p. This allows 
us to count modes that lie outside the horizon and below the momentum 
cutoff. 

The relevant mode density is the inverse of the frequency gap Au; that 
separates two neighboring modes which are orthogonal with respect to the 
one-particle Hilbert space inner product restricted to the black hole exterior. 
As shown in the Appendix, in the WKB approximation, this mode density 
is given by 

n(u,£) = {2ttkY 1 ln(fc max /fc min ), (8) 

where /c max is the cutoff wavevector A and k m \ n is the value of k u} i(p) when 
the trajectory crosses the horizon after having reflected from the angular 
momentum barrier (see Fig. [I]). In the near horizon region, the wavevector 
decreases exponentially as exp(— kt) from /c max to k m \ n along the trajectory, 
thereby allowing to write (|8|) as n(u>,£) = At/2tt. 

We now explain why it suffices to restrict attention to the near horizon 
region. At the turning points of the classical trajectories the dispersion 
relation ([7]) implies 

where p is the transverse momentum defined by 



l-vL.~(u/p) 2 , (9) 



p 2 = £(£+l)/r 2 . (10) 

Since (as it turns out) the entropy is dominated by the high angular mo- 
mentum modes with p>K, and the contribution of high frequencies to the 
entropy is exponentially suppressed by exp(— 2ttlu/k), we have (1— v 2 )t p <C 1. 
It therefore suffices to use the linear approximation ([3]), and to set r = Tu 
in the definition (|10p of the transverse momentum p. 

The dispersion relation (JJJ) at the horizon yields k m - m = (p 2 — lo 2 )/2u>. 
The mode density © in the presence of the cut-off is therefore given by 

n cut (u,£) = (27TK)- 1 In (^^f) (11) 

The condition that k m i n be less than k max = A constrains p to be less than 

p max = ^2ujK + uj 2 , (12) 

thereby eliminating the quadratic divergence in ([5]) that comes from the 
angular momentum sum. For transverse momenta p close to to the contribu- 
tion to the entropy is not well determined by our near horizon approximation 



S 




Figure 1: Trajectories of the particle and partner in the relativistic and 
dispersive cases. The density of states is always proportional to the PG 
time lapse At. In the relativistic case (left) this is the time the particle 
spends outside the horizon and below the hard cutoff. In the super-luminal 
case (center) it is simply the time the particle spends outside the horizon. 
In the sub-luminal case (right) it is the time between the horizon crossings 
of the partner and the particle. 
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since, as shown by Q, the turning point is not near to the horizon. How- 
ever, these momenta make a negligible contribution to the entropy. Thus 
we simply adopt the lower limit p m i n = to that is imposed by the condition 
fcmin > 0, which is obeyed by any mode that is initially outgoing. 

Since the atmosphere entropy is dominated by large angular momenta 
we can replace the sum over I in ([5|) by an integral over p 2 , 

J2(2£ + l)^r 2 H [ dp 2 . (13) 

Then inserting (jlip for the mode density into © yields 

Satmo, cut = Ax («A) x (27r)~ 4 / dx X S thcrra (x) , (14) 

Jo 

where A is the horizon area 4irrfj and x = 2ttlo/k. (The integral is equal 
3£(3) ~ 3.6.) This expression gives the entropy of the black hole atmosphere 
when regulated by a hard momentum cut-off imposed in the freely falling 
frame that is at rest at infinity, evaluated using the WKB approximation. 
We now comment on several aspects of this result. 

Frame dependence — Had we imposed the cut-off A in a general free-fall frame 
the result would have been the same as (fl4j) . except with k replaced by k v 
defined in ([3]) or equivalently ([1]). Thus the entropy actually depends on the 
free-fall frame in which the cut-off is defined. This is to be expected, since 
the cut-off prescription is not Lorentz invariant. 

Universal dependence on transverse momentum cutoff — Although the en- 
tropy depends on the free-fall frame adopted for the cut-off A, it scales 
universally with the maximum angular momentum as 

I 2 ~ At? (15) 

^max ^Fmaxi \ xt -V 

with p max = \/2k v A the value of p max (fl2j) evaluated at to = k v (where most 
of the entropy resides). Thus each unit of area (p m ax) -2 carries of order one 
bit of entropy. Put differently, the entropy scales as the number of angular 
momentum states below the cutoff. The same scaling is found in the case of 
a dispersive cutoff in the next section. 

Near-horizon approximation — As explained after (|10p . the near- horizon ap- 
proximation is justified for modes with p»K. We can now see that this is 
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satisfied for the dominant modes, since p max ~ (kA) 1//2 3> n. 

Anisotropy of momentum cutoff — Since p max is much smaller than the cut-off 
A, the entropy would be unchanged were the cut-off A imposed isotropically 
on the norm of the momentum rather than on the radial component alone. 

Mode density for dominant angular momenta — Since the entropy scales as 
the number of angular momentum states, the number of modes contribut- 
ing to the entropy at fixed angular momentum (£, m) is of order one for 
the dominating values of I. Indeed, evaluating the mode density (jlip for 
p parametrized as a fraction of the maximum value, j? m ax/A^, yields n cu t = 
(In N)/2ttk. Since thermal suppression by the Boltmann factor limits the 
effective frequency range to Alo ~ k/2tt, the number of contributing modes 
is ~ n cut Auj ~ (In N)/(2ir) 2 . Since most of the entropy comes from N close 
to unity, only a "fraction of a mode" contributes for each angular momen- 
tum. In this regime ra cu tAu; should be interpreted as the probability to find 
the corresponding particle in the atmosphere outside the horizon. 

Comparison with "brick wall" model — The entropy (|14p scales as kA, linearly 
with the momentum cutoff A. This scaling differs from that of the entropy 
regulated with a "brick wall" cutoff [7], i.e. a reflecting boundary condition 
at a proper distance e from the horizon, which scales quadratically as e~ 2 . 
The difference arises only because different frames have been used to measure 
the cutoff scales: we measured A on a time slice orthogonal to the free-fall 
worldlines, whereas e was measured in Ref. [7J along a static time slice, 
orthogonal to the Killing vector. On the free-fall time slice, the same brick 
wall lies at a much smaller proper distance 5 = ne 2 /2 from the horizon. 
Hence when expressed in terms of 5, the brick wall entropy scales as <5 _1 , in 
agreement with what we found. 

The agreement in overall scaling with the cutoff is a feature of the total 
entropy, after integrating over all frequencies and transverse momenta. This 
agreement arises because the dominant contributions scale the same way, 
despite the fact that the physical nature of the two cutoffs is different. In 
particular, the brick wall cuts off all modes at the same location, whereas 
we cut off all modes at the same momentum. In our case the cutoff location 
is not sharp, but can be characterized using the WKB approximation. Its 
value then depends on the frequency and angular momentum. The angular 
momentum dependence is negligible because, as noted after equation (fT2|) . 
p <C A for the frequencies u> ~ k that dominate the total entropy. The 
frequency dependence on the other hand is significant. Using the dispersion 
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relation ([7|) one finds near the horizon k ~ uj/ Vc(p — ph), so that the WKB 
position cutoff is uj/kA. However, this to dependence does not affect the 
total entropy scaling since frequencies much greater than k are thermally 
suppressed. A further negligible difference with the brick wall evaluation 
arises because we cut off the WKB trajectories at the final horizon crossing, 
whereas the brick wall trajectories begin and end at the same radius, as 
explained in more detail in the Appendix. 

4 Dispersion induced cut-off 

We now remove the high wavevector cutoff A and instead introduce disper- 
sion into the field propagation by including higher space derivative terms 
suppressed by powers of A. Both super-luminal and sub-luminal dispersion 
will be considered. As explained above, it did not matter much whether the 
hard cut-off was isotropic or only in the radial direction. In the dispersive 
case the same will be true, at least in the super-luminal case. Whether it 
is true in the sub-luminal case will depend upon an open question to be 
discussed below. 

The relativistic dispersion relation between frequency f2 and wavenumber 
q for a massless field in flat spacetime is Q 2 = q 2 , where q 2 = q • q. We 
consider an isotropic UV-modified dispersion relation of the form 

n 2 = F 2 (q) = (q±±q n /A n - 1 ) 2 , (16) 

chosen for convenience so that F 2 (q) is a perfect square. The quantity A 
sets the scale for the nonlinear dispersion, and we assume n > 1. At low 
wavenumbers q <C A the dispersion is relativistic. At high wavenumbers the 
group velocity v g = d q Q is super-luminal for the plus sign and sub-luminal 
for the minus sign. At q = A we have Q(A) = (1 ± ^)A and v g (A) = 1 ± 1. 
We now adopt units with A = 1. 

The dispersion relation has the given form only in a particular rest frame. 
As in the case of the hard cutoff we choose this to be the free-fall frame that 
is at rest at infinity. As in the nondispersive case, we can decompose the 
scalar field in the form ([6]) and use the WKB approximation for the radial 
functions u w f. The positive free fall frequency outgoing {k > 0) mode with 
angular momentum t satisfies the dispersion relation 

1 n 

n = u-v{p)k = (k 2 +p 2 )2 ±l(k 2 +p 2 )2 (17) 
« k + p 2 /2k±k n /n, (18) 
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where the approximate form in the last line holds when p <C k and k <C 
A = 1. As can be verified a posteriori, this approximation is valid in the 
present context, as in the hard cutoff case. Moreover, the r dependence in 
p 2 = 1(1 + 1)/V 2 can be safely replaced by the constant value rn- 

As before we take the field to be in the free fall vacuum state. A num- 
ber of studies have demonstrated that in this state, to leading order in k/A, 
models with UV dispersion reproduce the relativistic result that the outgoing 
modes at infinity are thermally populated at the Hawking temperature [301 
EU [321 [33] . To evaluate the entanglement entropy we need to know whether 
thermality extends to the near-horizon modes that contribute to the entropy. 
In the Appendix we deduce the quantum field theoretic description of the 
near-horizon free fall vacuum, generalizing to four dimensions the treatment 
of Ref. [34]. The state is a product of two- mode squeezed states that are 
thermal when reduced to one mode. We also discuss in the Appendix how 
the WKB approximation can be used to identify those entangled modes 
localized on either side of the horizon. In this section we shall assume that 
this WKB treatment provides an accurate evaluation of their density. 

Under the above assumptions, the atmosphere entropy takes the same 
form ([5]) as before, but with a different density of modes n^i Bp (u, £). As 
shown in the Appendix, the mode density is again given by ([8]). As in the 
relativistic case, the wavevector k m [ n corresponds to the point where the 
particle trajectory finally crosses the horizon. In the WKB approximation, 
the only place where the dispersion enters the calculation of the atmosphere 
entropy is in the determination the initial wavevector fe max . 

In the super-luminal case the outgoing positive frequency modes emerge 
from behind the horizon after parting ways with their negative frequency 
partners which always remain behind the horizon (see Fig. [1]) . Therefore the 
relevant trajectory for determining the density of states to be included in 
the entropy extends between the two horizon crossing points of the positive 
frequency trajectory, and fc max is the value of k at the first (outward) horizon 
crossing. The mode density at fixed u> and I is thus finite, and the sum 
over £ is cut off when the turning point drops below the horizon and the 
two horizon crossing points merge. The atmosphere entropy is thus finite. 
We now discuss the determination of k max in the sub-luminal case, before 
computing the entropy for both cases together. 

In the sub-luminal case, the positive and negative frequency outgoing 
modes are initially infalling with respect to the static frame. For high enough 
frequency or low enough angular momentum the positive frequency modes 
convert to outward propagation at a turning point just outside the horizon 
(see Fig. [I]), followed by a second turning point at the angular momentum 
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barrier where they begin to fall back in. For lower frequency or higher 
angular momentum the pair of turning points disappears and the particle 
just follows the partner straight across the horizon with a small time lag. 
Since the entanglement of these modes is the source of the entropy, we 
shall identify fc max the value of k for the positive frequency particle at the 
time the negative frequency partner crosses the horizon, leaving the particle 
stranded outside. This notion requires us to determine what is the particle 
momentum at a given PG time when the partner has a given momentum k. 
This determination is impossible if we refer only to Hamilton's equations, 
but in quantum mechanics the positive and negative frequency modes are 
linked through the specification of the in-modes which define the vacuum 
state. As shown in the appendix, by building wave packets of in-modes one 
can see that the particle and partner possess the same mean momentum at 
the same PG time. We therefore identify /c max as the value of the momentum 
when the partner crosses the horizon, see Fig. [TJ 

The mode density n(uj,£) of ([8]) is thus again finite, since A; max /&min 
is finite. There remains a potential divergence in the sum over angular 
momenta however since, unlike in the super-luminal case, this ratio remains 
strictly greater than unity. This is because the partner always crosses the 
horizon a finite time before the particle, as explained in the Appendix. The 
mode density therefore remains non-zero no matter how large i becomes, 
and it may not approach zero fast enough for the sum over £ to converge. 
We will examine this potential divergence more closely later on. 

The values of /c max and k m \ n are given by appropriate roots of the dis- 
persion relation (fl~8j) evaluated at the horizon, 

u = p 2 /2k ± k n /n. (19) 

The upper and lower signs are for the super- and sub-luminal cases respec- 
tively. This can be expressed in terms of the scaled variable k = k/\u>\™ 
as 

akT~k n+1 =p 2 /2|w| 1+ i (20) 

where a is the sign of co, positive for particle and negative for partner. The 
ratio fc max /fc m i n is the same as the ratio of the fc's, hence it depends upon u> 
and p only through the combination p 2 . 

Following the above discussion, the mode density n(u>,£) in the super- 
luminal case is determined by the log of the ratio of the two roots fc max 
and fcmin of ([20]) for a = +1 (particles) and the upper (negative) sign. The 
maximum transverse momentum is determined by the maximum of the left 
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hand side of ([20]) . It is given by 



where the explicit dependence on the momentum scale A has been restored. 
This is essentially the same as in the case of the hard cutoff (|12p in the limit 
n — > oo, when the dispersion turns on suddenly at k = A. When n is finite, 
Pmax scales with a lower power of the "new physics" scale A than for the 
hard cutoff. Therefore the atmosphere entropy will be similarly reduced. 

In the sub-luminal case the lower sign applies, and /c m ; n is the root for 
a = +1 (particle) while fc max is the root for a = — 1 (partner). In both cases, 
the mode density can be expressed as 

1 1 ( P 2 \ 

ndisp(^) = — ln(fc max /fc min ) = - — Af n —j , (22) 

where is a function of its dimensionless argument that cannot be written 
in closed form. In Fig. [2] we plot this function for both sub-luminal 
and super-luminal (M^) dispersion in the case n = 3. Notice that while the 
cutoff scale is taken to be A = 10 4 k, the maximal value of the transverse 
momentum in the super-luminal case is only ~ 25k. 

Using the mode density just obtained, the entropy (J3J) (again using (Tl3|) ) 
evaluates to 

/ K \ n I 

Satmo,di S p = A X KA X f - J xC±, (23) 

where the dependence on A has been restored. The factor is a constant 
defined by 

, 1 , rVmax,n,± , \ / poo , . 

Ct = {2n)- i ~-[J o dyN±(y))(j o dx x 1+ ~ s thcrm (x)) , (24) 

where y is the dimensionless variable p 2 /(2ojA(uj/A) n ). The upper limit 
of integration y ma x,n,+ in the super-luminal case is the maximum value of 
the left hand side of (|20p . yielding a finite total entropy within the WKB 
treatment we have adopted. In the sub-luminal case j/ maXjn _ = oo, as 
explained above, hence the result depends on the asymptotic mode density 
at momenta larger than A. 

The expression (|23p gives the atmosphere entropy in the dispersive case. 
As in the case of the hard cutoff, the entropy scales universally with the 
square of the angular momentum cutoff, i.e. as ^ max ~ ^4p max . In the 
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Figure 2: Plot of the mode density (|22|) vs. the transverse momentum p 
at fixed frequency u = K, for sub-luminal (upper curve) and super-luminal 
(lower curve) dispersion, with A/k = 10 4 and n = 3, in units with k = 1. 



present case, p max is the value of (|2ip at w = k (where most of the entropy 
resides). Comparing with our previous result (114p for the case of a hard 
cutoff, we note that the scale A enters with a different power which depends 
on n. In the limit n — > oo one recovers the hard cutoff result, at least in the 
super-luminal case. This happens because for large n the dispersion turns 
on abruptly at k = A. 

Asymptotic mode density in sub-luminal case — The absence of an upper 
limit to the angular momentum sum in the sub-luminal case means that the 
result depends on the behavior of the dispersion relation at momenta higher 
than A. The dispersion relation (|16p was introduced to parametrize only 
the first order correction to the relativistic dispersion. To analyze the mode 
density in the regime of very high momenta the first term in the Taylor 
expansion is inadequate, so we allow for a general dispersion relation of the 
form (u — vk) 2 = F 2 where F(k,p) is now an arbitrary function giving rise 
to sub-luminal propagation. 

At fixed k and p, as shown in the Appendix, the dispersion relation 
implies that the particle and partner locations near the horizon are sepa- 
rated by Ap = 2uj/(nk), with the particle at larger p. The mode density 
is determined by the time between the horizon crossings of the partner and 



16 



the particle, At |Ap/(d/9/dr)| = 2u/(Kk(l — F^))- It is hence given 
approximately by 



where k satisfies the dispersion relation at the horizon, uj + k = F(k,p). 
As long as F ^> u>, we have k — F(k,p) hence n su b behaves as ~ uj/F 
since for the sub-luminal case F^ is assumed less than one. The number 
of modes in the relevant thermal frequency range Atu ~ k is therefore ~ 
k/F <C 1, assuming F 3> k. That is, there is only a tiny probability to 
find a particle with large transverse momentum. The WKB approximation 
is probably unreliable in this regime, so our methods are unable to evaluate 
the contribution of these modes to the entropy. 

5 Discussion 

We have found that a freely falling hard momentum cutoff produces an 
entanglement entropy that scales as A k v A with the cutoff A and the surface 
gravity k v defined in a given free fall frame. When the frame is taken to be 
at rest at infinity k v coincides with the usual surface gravity k. With no 
cutoff but instead sub- or super-luminal dispersion with index n the scaling 
has an additional factor of (ft/A) 1 /™. If the freely falling cutoff A is of 
the order of the Planck mass M p this entropy is negligible when compared 
to the Bekenstein-Hawking entropy Sbu ~ AMp 2 (for macroscopic black 
holes). This would conveniently eliminate the species problem, i.e. the non- 
universal contribution of the entanglement entropy to the black hole entropy, 
without the need to suppose a relation between the entanglement entropy 
and renormalization of the low energy effective Newton constant. 

If instead the cutoff were given by A = Mp /k in the freely falling frame 
the scaling of the entanglement entropy would agree with <Sbh- This is 
equivalent to saying that the cutoff is actually M p in the static frame. While 
this might be considered natural from a "geometric" viewpoint, it seems 
unnatural from the viewpoint of the equivalence principle, which suggests 
that the cutoff should be a property of the freely falling vacuum. However, 
the equivalence principle would not forbid a cutoff arising from quantum 
gravity effects sensitive to the global properties of the background. In fact, 
there is some evidence for this. 

Strong gravitational interactions between ingoing and outgoing quanta 
arise when the invariant mass squared reaches Mp [7] . The radiative correc- 
tions associated with these interactions in the free fall vacuum outside the 




(25) 
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black hole were studied in Ref. \17\ [T8] . Beyond the local effects that would 
also be present in the vacuum of flat spacetime, it was found that there 
are further contributions, dominated by in-modes with Killing frequency 
Wi n ~ «■ These corrections induce a dissipating cutoff on the (backward) 
propagation of outgoing modes at the location where their frequency is / k 
in the freely falling frame, i.e. M p in the static frame. While the significance 
of that preliminary analysis remains unclear, it does indicate a way in which 
the cutoff entanglement entropy might account for the Bekenstein-Hawking 
entropy. 

It would be interesting to extend the considerations of the present paper 
to the cases of de Sitter and Rindler horizons. In the near- horizon limit these 
are identical to a black hole horizon, and a strong case can be made that 
they too posses a universal entropy density l/4Lp [35J. On the other hand, 
they differ in the nature of their "observer-dependence" and, perhaps more 
importantly, in the nature of the globally stationary free-fall frames that 
they admit. We leave the study of the regularized entanglement entropy for 
such horizons to future work. 

We would like to reiterate that the WKB approach employed in this 
paper cannot provide conclusive results regarding the finiteness or scaling 
of the entanglement entropy in the dispersive cases, and it is particularly 
inconclusive in the sub-luminal case. The uncertainty arises since the mode 
density for the relevant modes is smaller than unity, so that we are not in 
the regime where a WKB treatment is guaranteed to be reliable. Therefore 
only a field theory analysis of the reduced density matrix outside the horizon 
would permit definitive conclusions to be reached. 

In this paper we have studied the entanglement entropy between the 
degrees of freedom on either side of the sharp, low frequency horizon on a 
fixed black hole background. In the relativistic case, this horizon is a strict 
causal barrier, but in the dispersive cases it is dynamically distinguished 
only in the low frequency limit. (For sub-luminal dispersion the location of 
the barrier is "frequency dependent" , while in the super-luminal case there 
is no barrier at all.) While this entanglement entropy is interesting, a more 
important, physical problem is to evaluate the "actual" thermodynamic or 
statistical entropy. The physically relevant definition of entanglement is not 
so clear in the dispersive cases, and in fact is unclear even in the relativistic 
case, once quantum fluctuations of the horizon are admitted. Presumably, 
the definition should ultimately be provided by a specified set of physical 
observers, interacting with the system in a definite manner. 

Note added: While this paper was in press we became aware of a paper 
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by Frolov and Novikov [37] in which a regulated black hole entanglement 
entropy for a massless scalar field was computed. A method similar to 
the one used here using WKB trajectories was employed to identify the 
modes contributing to the entropy. The cutoff was imposed by including 
only contributions from modes outside a certain minimum distance from 
the horizon, motivated by computations of quantum horizon fluctuations. 
The minimum distance was taken to be of the order the Planck length as 
measured on a static slice, so that the entropy scales with the area in Planck 
units. 
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Appendix 

In this appendix we spell out the field theoretic rationale for our WKB evalu- 
ation of the mode density relevant to the entanglement entropy. The central 
ingredient is a description of the correlated structure of the near-horizon 
free-fall vacuum in terms of modes written in momentum space. These re- 
sults could be useful in a future approach to computing the entanglement 
entropy without recourse to the WKB approximation. 

Action & field equation 

We consider a massless scalar field in the spacetime with line element ([2]). 
It is governed by the action 




{d^f}- (A.l) 
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The field equation arising from this action is 

{d T + dpv)r 2 (d T +vdp)(p-d p r 2 dpLp--?— d e sm9d S (p- . \ - dip = 0. (A. 2) 

sin sin a r 

Modes 

Modes of fixed frequency and angular momentum take the form 

Vulm = exp(-iwr)r" 1 n^(/))Y £m (6',0). (A. 3) 
To find the equation satisfied by the mode function u^i we compute 

d p r 2 d p r- 1 = rd 2 p - r", (A.4) 

and 

(— iui + dpv)r 2 (— iu> + v<9 p )r _1 = r(—iu> + d p v)(—iu> + vd p ) — (v r')', (A. 5) 
where the prime ' denotes d p . The mode equation is thus 

{-iuo + d p v)(-iu + vd p ) -d 2 + Vt{p)}u u = 0, (A.6) 

with 

y<M = (i-V_W^ + <C£+i), (A . 7) 

This is written for the general line element ([2]). With the Painleve-Gullstrand 
(PG) coordinates for the Schwarzschild metric we would have r(p) = p, so 
the first term in Ve vanishes and the second term becomes just —2vv'/r. 

Near-horizon approximation 

We are interested in modes propagating very near the horizon since these 
are the ones that contribute most of the entanglement entropy. Expanding 
about the horizon in x = p — pu we have 

v(p) = -1 + k v x + --- (A. 8) 

where k v = v'(0) is the surface gravity. The higher order terms in x are 
negligible in the near- horizon region, and the 0{x) term is only significant 
when multiplied by a large coefficient. The derivatives d p in the first term 
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of (|A.6[) provide large coefficients for the modes of interest, due to the large 
redshift at the horizon, so we retain the 0(x) contribution to the velocity 
(|A.8|) in that term. Since the entanglement entropy is dominated by large 
angular momenta, we can drop from Vg all but the ^-dependent term. In the 
near-horizon region we can set this term equal to its value at the horizon, 
Vi = p 2 , where we have introduced the notation 

p 2 = e(£ + l)/r 2 H (A.10) 

for the transverse momentum squared evaluated at the horizon. 

Momentum space 

The momentum space modes u{k) are related to u(p) by Fourier transform, 

u(p) = (2vr)~ 1/2 J dk u(k) exp(ikp). (A.ll) 

In momentum space the mode equation ()A,6|) becomes 

(iv - kv)(uj - vk)u u i = (k 2 + Vi)u w i, (A.12) 

where v = v(p) and Vg = Vi(p) with p = idk- In the near-horizon approxi- 
mation, and for large i this mode equation takes the simpler form 

(oj - kv)(u - vfyu^i = F^tyuui, (A.13) 

where now 

v = — 1 + K v p = — 1 + iK v d p , (A. 14) 

and 

F 2 (k) = k 2 +p 2 . (A.15) 

To simplify the notation we now drop the subscript v on k. 
Thanks to the identities 

{uj-vk)k- iu}/K ~ l e- ik/K u = k~ lu) l K - l e- ik/K {-iKkd k u), (A.16) 
(w - kv)(to - vk)k-™' K - l e- ik ' K u = fc-*"/«-i e -«/" (- K 2 k 2 d 2 k u),(A.17) 

the w-dependence of the solutions to (|A.13h can be factored out, so the 
solutions take the form 

u ut = k~ iu /*- 1 e- ik / K xt(k), (A.18) 
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where xi obeys 



-K 2 k 2 d 2 kX t = Fixe, (A.19) 
= (k 2 +p 2 )Xi- 

For k 3> p the solutions are exp(±ik/K), while for k <C p they are k a with 
a = | ± \J\~ P 2 /k 2 - Only the first case is relevant for the modes outside 
the horizon however, since k m - m ~ p 2 /<-o which is much greater than p for the 
modes relevant to the entropy. 

It so happens that (|A.19P is closely related to Bessel's equation, and is 
exactly solvable in terms of Bessel functions: 

Xi(k) = VkC v {k), v = ^\-p 2 Ik 2 (A.20) 

where C = J,Y, H^' 2 \ As shown below these solutions can be classified 
according to whether xe has positive or negative Wronskian 

W = -i(Xid k X£-(dkXi)Xe)- (A.21) 

Dispersion 

Now we consider modifying the field equation by the addition of higher 
spatial derivative terms in the local rest frame defined by the unit timelike 
vector field w = d T + v (p)d p . For the metric ([2]) this frame is orthogonal to 
the surfaces of constant r. 

The dispersion of interest becomes important at high wavevectors, where 
"high" is specified by a momentum scale A>k, Since the modes of interest 
are only important near the horizon, it is sufficient to study dispersion by 
a modification of the near-horizon mode equation. To be able to handle 
arbitrary dispersion, it is convenient to work in momentum space |31| 134] . 
Allowing for distinct modifications of (|A. 13|) in the radial and transverse 
directions yields the modified near-horizon mode equation 

(uj - kv){uj - vk)uu P = F 2 (k,p; A)n wp , (A. 22) 

where F(k,p;A) is a function determining the dispersion, which is assumed 
to be relativistic if k,p <C A. 

As before, the w-dependence of the solutions to (|A.22p can be totally 
factored out, so the solutions again take the form (IA.18jl . 

Uu P = k-^/^e-^^XvAk), (A.23) 
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but with xe replaced by Xp.A- The latter satisfies (|A.19|) with the new F 2 
function: 

- K 2 k 2 d 2 k x P ,A = F 2 (k;p,A) XpA . (A.24) 
Klein-Gordon inner product 

On a constant r slice the KG inner product in both the relativistic case and 
in the presence of dispersion is given by 

{vl,P2)=i J dttdpr 2 (ip*(d T + vd p )(p 2 - {{d T + vd p )tp*)(p2), (A.25) 

because the dispersion we consider only affects the spatial part of the action. 
It is diagonal on angular momentum modes, and for partial waves © with 
the same angular momentum it takes the form 

(<Pu>',fu) = J dp(uu,(u + ivd p )u u + ((u' - ivd p )u^)*u^. (A.26) 

In momentum space it becomes 

(<Pv',<Pu) = J dk(ul>(u> - vk)u u + ((a;' - vk)^)* . (A.27) 

For solutions factorized as in (|A.18j) or in (|A.23|) . we make use of (|A.16|) to 
obtain 

((Put,?*) = -»« / dkV* u / K - 1 ( X *d kX -(dkX*)x) (A-28) 

where Aoj = oj 1 — u. In addition, for any solutions of (|A.19|) or its dispersive 
modification (|A.24|) . the Wronskian of X is a constant. Thus, irrespective of 
angular momentum, dispersion, or WKB approximation, the KG overlap of 
two modes in the near-horizon region is proportional to 

J dkk^l K - x = J d(\nk)e^ lnk ^ Auj / K . (A.29) 
Mode density and classical trajectories 

To compute the entanglement entropy we need to know the density n(uj,£) 
of modes exterior to the horizon. This is just the reciprocal of the gap 
in Killing frequency u that separates one exterior mode from the nearest 
orthogonal exterior mode with the same angular momentum. Orthogonality 
here is defined with respect to the one-particle Hilbert space inner product, 
which is defined by the Klein-Gordon inner product. The mode density is 
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thus the inverse of the minimal Aw for which the overlap (|A.29P vanishes, 
hence 

n(u, i) = (27TK)- 1 ln(fc max /£; min ). (A.30) 

We implement the restriction to exterior modes through the limits of inte- 
gration /c max and fc m in, which are determined using the WKB approximation. 

In both the relativistic and dispersive cases k m { n is the value of the 
WKB momentum k u i{pn) at the horizon crossing where the trajectory of 
the particle falls back into the black hole after having reflected from the 
angular momentum barrier (see Fig. [I]). The value of & max is just the hard 
cutoff A in the relativistic case, and in the dispersive case it is determined by 
the passage of the particle (or its partner) across the horizon, as explained 
in the text. 

The dispersion relation associated with the mode equation (|A.22p is 

to — v(p)k = F(k,p; A), (A.31) 

where we took the positive root to describe outgoing modes. The classical 
trajectories satisfy Hamilton's equations, 

dp/dr = du/dk, (A.32) 
dk/dr = -du/dp, (A.33) 

where the Hamiltonian is given by u(p, k; £), solution of (|A.3ip , In the near- 
horizon region, where v (p) ~ — 1 + kx, (|A.33|) implies the simple relation 

dk/dr = -nk, (A.34) 

both in the relativistic and dispersive cases. The wavevector k therefore 
decreases exponentially as exp(— hit) along the trajectory, both before and 
after reflection from the angular momentum barrier (which for large I lies 
in the near-horizon region). Notice that the same exponential law with 
k > applies to the trajectories of the (outgoing but trapped) partners 
characterized by the opposite Killing frequency — u. The sign of k should 
therefore be thought as that of the FF frequency. 

Using (|A.34|) . the mode density (|A.30P can thus be reexpressed in the 
simple form 

n(w,f) = (2^ 1 Ar(w,f), (A.35) 

where At(lo,£) is the lapse of PG time evaluated along the relevant part of 
the classical trajectory. Two near-horizon modes are therefore orthogonal 
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with respect to the KG product if their frequencies are separated by Aw = 
2it/At, irrespectively of the value of p and the dispersion law -Fill 

It is informative to compare our result with the mode density for the 
static "brick wall" reflecting boundary condition used by 't Hooft in Ref. [7]. 
In that model, the quantum field is required to vanish at the wall, i.e. 
at a fixed radius outside the horizon, 't Hooft expressed the mode den- 
sity as the derivative of the number of radial nodes with respect to fre- 
quency, vr -1 / dr (dk r /duj), where to is the frequency with respect to the 
Schwarzschild time coordinate. Using the group velocity relation dr/dt = 
diojdk r this becomes simply 7r _1 J dt. The integral is the Schwarzschild time 
lapse along the classical trajectory that rises from the wall and stops at a 
turning point. 

The one-bounce trajectory that rises and then falls back to the wall is 
time-symmetric in terms of the Schwarzschild time coordinate, so we may 
extend the integral over t to the full round trip integral and divide by two. 
This results in an expression that looks identical to (|A.35p . but it differs in 
that the time coordinate is the Schwarzschild time rather than the PG r 
coordinate in ([2]). However, if we choose v(oo) = 0, then these time 

coordinates differ only by a function of radius, so their lapses at any fixed 
radius are equal. A small discrepancy remains because, unlike t 'Hooft 's 
brick wall, our regularization scheme leads to a slightly asymmetric domain 
of integration. The final point of the trajectory defining our mode density is 
at horizon crossing on the inward branch, while the initial point lies outside 
the horizon where the momentum is equal to the cutoff. The lapse of time 
associated with the extra bit is of order 1/A however, and is thus negligi- 
ble for our purposes. Our cutoff mode density and that of the brick wall 
prescription would thus agree if the cutoff momentum occurred at the same 
location as the brick wall. The comparison of the brick wall and momentum 
cutoffs is discussed further in the text. 

7 The linear growth in time of the relevant mode density in (| A.35[) is familiar from 
the behavior of E mo( j es = (Ai/27r) J du>p(uj), the number of modes involved in transition 
probabilities (the Golden Rule). Less familiar is the peculiar fact that here this appears 
expressed as (|A.30|) . coming from an integral over momenta of the form J dk/k. This 
feature of near horizon states was previously noticed in [36] in the context of both the 
acceleration horizon of an Unruh detector, and when counting the number of Hawking 
quanta in terms of the free fall frequency they had near the horizon. See the discussion 
around Eq. (3.44) in [25] for a brief account of this. 
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Free-fall vacuum: localization of particles and partners 



It follows from (|A.28|) that the modes 

= 0(k)k~ iu ' K - 1 e- ik '"x£(k), (A.36) 

have positive norm if xi nas a positive Wronskian. These are complete 
and orthogonal when u ranges over all values, both positive and negative. 
They thus generate the Hilbert space of one-particle outgoing states. The 
outgoing free-fall vacuum is defined by the condition that it is annihilated 
by the destruction operator = (4>^£,4>) corresponding to these modes 
because their FF frequency f2 = u — vk is positive, see (|A,16|) . 

The Fourier transforms of the modes (|A.36p have support on both sides 
of the horizon, even for £ = 0. At the quantum level this gives rise to 
correlations between field fluctuations on either side of the horizon in the 
free-fall vacuum. These correlations are best characterized using a different 
mode basis consisting of positive Killing frequency "particles" and negative 
Killing frequency "partners" , as we now explain. 

The key mathematical ingredient is the decomposition into terms ana- 
lytic in the upper and lower half complex £>planes: 

9(k)k~ ilJJ ~ 1 = 1 [e w (lfe - ie)~ iw - 1 - e~™(k + ie)""^ 1 ] (A.37) 

where e — * + . Here we have chosen the branch cut of the logarithm to run 
along the negative k axis. (To simplify the notation we now adopt units with 
K = 1.) The Fourier transform of the first term vanishes for x < 0, since 
it is analytic in the lower half /c-plane and bounded as \k\ —* oo. Similarly 
the Fourier transform of the second term vanishes for x > 0. This is not the 
whole story however, since the mode (|A.36|) contains the additional factor 
e- ik Xi(k)B 

Let us first consider the localization in the relativistic case with vanishing 
angular momentum. In this case X£=o(k) oc e ±ifc . The solution with positive 
Wronskian is e tk , which cancels the factor e~ tk in the mode. Thus the mode 

8 For zero angular momentum modes xo = exp(±ifc/fc) the sign of the Wronskian is ±, 
and (|A.16|) becomes (u) — vk)u^o — ifcit^o. Those modes with positive free-fall frequency 
and positive Wronskian therefore have only positive wave-vectors, i.e. they are outgoing. 
In the presence of small (large) angular momentum the would-be purely outgoing modes 
are partly (totally) reflected by the angular momentum barrier and fall back as ingoing 
modes. We can still identify these as "outgoing" however, in the sense that they have 
no incoming part far from the horizon. These are the modes that are relevant to the 
entanglement entropy. Those with large I constitute the dominant ones. 
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splits into a sum of two functions strictly localized inside and outside the 
horizon. If we allow for angular momentum, then the solution Xl(k) with 
positive Wronskian is no longer just e lk , but is rather given by the Hankel 
function, k^^H^ (h) in (|A.20p . This introduces branch cuts which lead to 
the Fourier transform of (k — ze) _,,w ~ 1 e~ lk k 1 / 2 Hu \k) having support inside 
as well as outside the horizon. The part inside corresponds to the reflected 
wave (of low momentum k ~ k) that has bounced off the angular momentum 
barrier. 

The addition of dispersion at a scale A further complicates the localiza- 
tion properties by modifying this time the UV part k ~ A S> n of the waves. 
As shown in |25} [M] the story is qualitatively similar, except for a mixing 
of outgoing and ingoing modes that corresponds to an extra pair creation 
process arising from the UV propagation of the modes. The amplitude of 
this is exponentially suppressed by a factor ~ e _7rA / K in the case of quartic 
dispersion F 2 = k 2 + k 4 / A 2 , and can therefore be safely neglected for black 
holes with surface gravity much less than A. 

When (|A.37j) is inserted into the mode expression (EOHl) the result ine 
two terms have KG norms with opposite signs. This is best seen with the 
help of the identity 

(k ± ie)-*"" 1 = O^k-^- 1 - 8{-k)e ±7TU) {-k)- iu >-\ (A.38) 

The contribution from 9(k) has positive norm while that from 6(—k) has 
negative norm, opposite aside from the factor e ±7rw . Hence (k — ie)~ l ^~ l 
yields a positive norm mode if oj > 0, while {k + ie)~ tU3 ~ l yields a negative 
norm mode. Taking into account the normalization factor of each of these 
modes, the annihilation operator for (|A.36|) with uj > thus decomposes as 

Vl _\-2,J ^-^oU- (A.39) 

That is, it is a combination of the annihilator of the (k — ie)~ lu '~ 1 mode and 
the creation operator for the complex conjugate of the (k + ie)~ lu) ~ l mode, 
which is a positive norm mode with Killing frequency -co. When to < 0, 
one finds (|A.39|) with b and c interchanged, and u replaced by —to. The 
condition that the annihilator for (|A.36p (for all to) annihilates the state thus 
becomes a statement that the state contains correlated pairs of b, c particles 
with opposite Killing frequency. Tracing over the negative Killing frequency 
particles yields a thermal state at the Hawking temperature Tu = k/2tt. 
Thus, even in the dispersive case, the near horizon free fall vacuum has a 
thermal character. 
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With vanishing angular momentum and no dispersion, these correlated 
pairs are strictly localized, with the positive Killing frequency one outside 
and the negative Killing frequency one inside the horizon. With angular 
momentum added, the positive Killing frequency mode falls back across the 
horizon with a low radial momentum. With dispersion, there is a further 
blurring of the outside/inside distinction which this time comes from the UV 
domain. Nevertheless, because A 3> k there exists an intermediate domain 
A 3> k S> k wherein (|A.37|) can be used to define a decomposition into a 
correlated pair of positive and negative Killing frequency modes which are 
localized on either side of the horizon. Moreover, when wavepackets are 
formed by integrating over u, those wavepackets will spend some time near 
the horizon behaving as do the ones with no angular momentum and no 
dispersion. The amount of time the one with positive to will spend outside 
the horizon while its partner is inside is (by definition) given by (|A.35|) . 

Indeed, in identifying the modes contributing to entanglement entropy, 
we use the WKB description of these modes, and require that the positive 
Killing frequency particle lies outside the horizon while the negative Killing 
frequency partner lies inside. Having formed a wavepacket, at any given 
time the wavevector of the particle and partner is given by the stationary 
point of the wavepacket integral over uj. It follows from (|A.37jl and the fact 
that no other rapidly oscillating function of to appears in the wavepacket 
that, at leading WKB order, this wavevector is the same for the particle 
and partner. That common wavevector is then used to locate the particle 
and partner using the dispersion relation with to > for the particle and 
— lo for the partner. 

Particle and partner trajectories 

In the near horizon region, the dispersion relation (1A.31|) reads 

u = Kkx - k + F(k,p;A). (A.40) 
Solving for the trajectory x = x(k,iv,p; A) we have 

x(to,p, k; A) = xo(k,p; A) + uj/ nk, (A. 41) 

with 

x (k,p;A) = (k-F(k,p;A))/Kk. (A.42) 

As just explained, the particle and partner have the same value of k at any 
given time, and opposite values of u. Hence the mean position (|A.42p of the 
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pair is independent of to, while the separation of the particle and partner is 
universally given by 

Ax = 2to/Kk, (A.43) 

independent of p and the form of the dispersion. This reflects at the classical 
level the fact, seen in ()A.23p . that the modes always split into a factor 
independent of p and F and a factor independent of to. This result, as many 
others, is identical to what was found in the 1+1 case studied in [33] because 
the transverse momentum p only enters as a modification of the dispersion 
function F. 

The space time trajectories x(t) are determined by (|A.41|) since k{r) is 
known. As seen in (IA.34[) . the fc-trajectories are simply given by k(r) = 
&oexp(— kt), again independently of the value of p and the dispersion rela- 
tion governed by F. 
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